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Abstract. We present new results on the decay patterns of scalar and tensor glueballs in the top-down holographic Witten- 
Sakai-Sugimoto model. This model, which has only one free dimensionless parameter, gives semi-quantitative predictions 
for the vector meson spectrum, their decay widths, and also a gluon condensate in agreement with SVZ sum rules. The 
holographic predictions for scalar glueball decay rates are compared with experimental data for the widely discussed gluon 
candidates /o(1500) and /o(1710). 
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INTRODUCTION 

Glueballs, color-neutral bound states of gluons, are the only physical states available in pure Yang-Mills theory and are 
also expected to show up in the meson spectrum of quantum chromodynamics (QCD) [1, 2, 3]. The mass of the lowest 
glueball state with - O'*'^ is predicted by lattice gauge theory [4, 5, 6, 7] to be in the range 1.5-1.8 GeV. However, 
as reviewed in Refs. [8, 9, 10], the identification of glueballs in the meson spectrum has remained elusive and will be 
in the focus of the PANDA experiment at FAIR [11, 12]. Various scenarios have been proposed, from unobservably 
broad glueball resonances [13, 14] to narrow glueballs that are identified with a dominant component of one of the 
fairly narrow isoscalar mesons in the mass range of the lattice results, /o(1500) or /o(1710) [15, 16, 17, 18, 19, 20, 21]. 

QCD sum rules seem to require a broad scalar glueball [13, 21], whereas the usually quite reliable large-Ac limit 
would indicate narrow states. The situation is similarly difficult in the case of tensor glueballs, which lattice simulations 
predict to occur in the range 2.2-2.6 GeV, but a clear identification of a corresponding /2 meson is missing. 

More information from first-principle approaches on the properties of glueballs would clearly be helpful. Lattice 
gauge theory should also be able to provide information on decay rates, but the extraction of real-time quantities is 
difficult, and glueballs are particularly difficult to investigate in the presence of dynamical quarks. 


HOLOGRAPHIC GLUEBALL SPECTRUM IN THE WITTEN MODEL 

A new approach to study glueballs and their properties has become available with the advent of gauge-gravity 
constructions where the supersymmetry and conformal symmetry of the underlying AdS/CFT correspondence is 
broken. The first and still most important example of a top-down holographic approach towards (large-Ac) QCD is 
the Witten model [22, 23] which is based on Nc D4 branes in type-IIA supergravity where one spatial dimension is 
compactified on a Kaluza-Klein circle with antiperiodic boundary conditions for the fermionic gauginos. This not only 
renders the gauginos massive, but through loop corrections also all adjoint scalar matter, so that in the limit where the 
Kaluza-Klein mass is sent to infinity one would arrive at pure-glue large-Ac Yang-Mills theory. Unfortunately, this 
limit is not accessible without leaving the realm of supergravity and going over to a full gauge-string duality which is 
too difficult for being of practical use. However, the hope is that already the supergravity approximation allows one to 
study nonperturbative features of large-Ac gauge theories qualitatively and that perhaps certain quantitative results are 
not too far from those of actual QCD, which is indeed frequently approximated quite well by large-Ac results. 

The supergravity background of the Witten model can be obtained by another dimensional reduction, namely a 
supersymmetry-preserving circle compactification of 11-dimensional supergravity with x" ^ x" -i-27rRii, which turns 
the AdS/CFT correspondence available for M5 branes to a nonconformal gauge-gravity correspondence for D4 branes. 
After introducing the supersymmetry-breaking circle compactification x‘* ^ x^ H- 27rR4, R 4 = the gravitational 
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FIGURE 1. The holographic glueball spectrum of the Witten model (left panel, masses in units of Mkk> “exotic” scalar mode 
in green) compared to the recent large-A^c lattice results of Ref. [24] (right panel, masses in units of the square root of the string 
tension). Dotted lines give the holographic glueball masses divided by from (11) with A = 16.63. 


background is given by a doubly Wick-rotated AdS? black hole geometry with (11-dimensional) metric 

ds^ = -jj[f{r)dxl + T]^ydx^dx''+ dx\^+—j^ + —dD}^, 

f(r)-l —rKK=^ATKK, (1) 

r° 3 

where L is the curvature radius of the asymptotic AdS? space forming a product space with an 5^ of radius L/2. With 
x^ and v" compactified, the boundary theory at sufficiently low energies is 4-dimensional Yang-Mills theory with 
parameters 

fl2 

^ym = = IngJsMKK, (L/2f = (2) 

Here gyM is the Yang-Mills coupling at the scale Mkk, normalized such that the QCD coupling as = g^ j{An) - 
g^y^l{2n) (see Ref. [25] for a discussion of this point). 

One of the earliest applications of gauge/gravity duality has been the calculation of the glueball spectrum from the 
supergravity background of the Witten model [26, 27, 28, 29, 30]. For this one needs to determine the spectrum of 
fluctuations of the metric and the 3-form field of the 11-dimensional supergravity theory. When classified according to 
their symmetries in the noncompact dimensions x^,...,x^ of the boundary theory this yields towers of scalar, vector, 
and tensor glueballs. This has been worked out completely in Ref. [30], with results displayed in Fig. 1 and compared 
with recent lattice results for the glueball spectrum in large-Ac gauge theory. Simpler (bottom-up) models [31, 32, 33] 
typically have fewer towers of glueballs, which include scalar glueballs dual to the dilaton, and tensor glueballs dual 
to the metric field. 

In the Witten model, the dilaton mode of type-IIA supergravity arises as 

3L^ 

dO=^dGii,ii. (3) 
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A scalar mode is obtained by having fluctuations of the form 





HD(r)GD(x), 


L2 


Hoir) 



df,dy 

□ 


Gd(x), 


(4) 


which also involves the trace part of metric fluctuations. The mass of this mode turns out to be degenerate with tensor 
modes from transverse-traceless fluctuations of only SG^y. 

As noted first in Ref. [29], the Witten model has another scalar mode associated with fluctuations involving also 
SGu, a metric fluctuation with components pertaining to the compactification dimension; 


6 Gaa = -j^fHG(r)G(x) 


SG^y = 


3r^ 

-^'kk 


4 5^6 _2^6 


df,dy 


KK/ 


Ml 


G(x), 


r 2 I 

L 
L2 

6 Grr^-^f 


SGn,n = j^-HG(r)G(x), 


-1 


3r^ 

-^'kk 


5r6-24^ 


HG(r)G(x), 


9Qr'^rl^ 

~ M|L2(5^6_2^6^)2 


HG(r)dM(x). 


(5) 


Ref. [30] subsequently found that this “exotic polarization” [29] has the smallest mass of all the supergravity modes. 

As displayed in Fig. 1, this does make the overall spectrum of the holographic glueballs appear to qualitatively agree 
with what is found in lattice gauge theory, but it leaves a certain abundance of scalar modes in the = O^'*' sector. 
Because the mass scale Mkk is not fixed in the Witten model, a more quantitative comparison with the lattice results 
is not possible. However, this changes when one extends the model to incorporate also quarks. 


INTERACTIONS WITH CHIRAL QUARKS AND MESONS IN THE 
SAKAI-SUGIMOTO MODEL 

In 2004, the Witten model was extended by Sakai and Sugimoto [34, 35] by a configuration of probe D8 and anti-D8 
branes that fill all of the 10-dimensional space of the Witten model except for the circle along x^. This introduces 
Nf «c Nc chiral quarks and antiquarks that are localized in x^, but since the D8 and anti-D8 branes have nowhere to 
end in the subspace spanned by and the radial (holographic) direction, these branes have to connect in the bulk 
geometry, leading to a simple geometric realization of chiral symmetry breaking (see Fig. 2). 

The low-energy theory of the (Witten-)Sakai-Sugimoto model involves a chiral Lagrangian for the Nambu- 
Goldstone bosons with a Skyrme term, which reflects the presence of a baryon vertex provided by additional D4 
branes wrapping the 8 “^ with Nc units of Ramond flux. It also provides the correct Wess-Zumino-Witten term re¬ 
quired by the chiral anomaly of QCD. With a minimal set of parameters, it represents arguably the most intriguing 
model of holographic QCD. As reviewed in Ref. [25], several features of the low-energy spectrum of QCD are nicely 
reproduced. 

Here we shall discuss in particular the predictions of the Witten-Sakai-Sugimoto model for the interactions of 
glueballs with mesons composed of chiral quarks, which were first studied in [36] and revisited by us in Ref. [37, 38]. 

Such interactions are determined by the action of the flavor D8 and anti-D8 branes, 

5d8 = -^DsTr J"d^xe~‘^ -\j- det (gMN + '^na'pMN) + Scs 

= -(2na'fTD^Tr jd\e-^^P§{l + ^~g^‘^~gQ^FpQFRs+0(F^)yScs ( 6 ) 
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FIGURE 2. Left: Configuration of the stacks of Nc D4 branes and Nf D8 and D8 branes in the 10 dimensions xo,...,xg, with 
X 4 taken as periodic. Right: Geometry near the horizon of the stack of Nc D4 branes, Nc » Nf, with u the radial coordinate in the 
transverse space JC 5 . 9 . [Figure courtesy of Andreas Schmitt] 


where the leading-order terms are explicitly given by 


c(F^) _ 
*^08 “ 



+ Ml^jf^FpzFyz 


^ ANc 
2167r^ ’ 


(7) 


with K=l + 7 }, where Z is a dimensionless holographic coordinate that runs from -00 to 00 as the connected branes 
(antipodal in x^) are followed from the boundary to the minimal value of the holographic radial coordinate r and back 
along the oppositely charged brane to the boundary. 

The global chiral symmetry S\J(Nf)xS\J(Nf) corresponds to a gauge symmetry with flavor gauge fields on the D8 
branes. Its spontaneous breaking is reflected by a massless mode Az - (pQ{Z)n{x) containing the Nambu-Goldstone 
bosons; massive vector and axial vector mesons are described by the even and odd eigenmodes of = iJ/„{Z)v^p\x). 
The lowest eigenvalue gives the identification 

=0.699314... (8) 


The real-world value of ntp can therefore be used to set the overall mass scale of the Witten-Sakai-Sugimoto model, 
leading to Mkk = 949 MeV. The masses of higher vector and axial vector mesons as predicted by the eigenvalue 
equation for are in fact surprisingly close to experimental results. For example, the predicted ratios iria^ jiUp and 
OTp* jiUp deviate only by 4% and 10% from the measured values. This agreement may however be somewhat fortuitous, 
for recent lattice simulations [39] at large Nc, where the holographic model should actually become more reliable, lead 
to deviations by 21% and 16%, respectively. Still, this can be viewed as a remarkable success given that the Witten- 
Sakai-Sugimoto model should show deviations from large-A^ QCD at scales above M^k and that already Vp is above 
Mkk- 


Other predictions of the Witten-Sakai-Sugimoto model for Nc-3 QCD depend on the value of the ’t Hooft coupling 
at the scale Mkk above which the model turns into 5-dimensional super-Yang-Mills theory. In the original papers 
[34, 35] the pion decay constant 


r 

Jn 


SAjt^ 


ANcM, 


KK 


(9) 


was fixed to f„ 92.4MeV, which together with Mkk = 949 MeV determines the constant k in (7) as 


K = dVJ(216;r^) = 7.45-10“^ 


( 10 ) 


corresponding to d = 16.63. 

As reviewed in Ref. [25], most other possibilities of fixing A point to somewhat smaller values. As an alternative we 
shall consider the recent large-Vc lattice result for the ratio Wp/ -\/cr = 1.504(50), where cr is the string tension, together 
with the holographic value 


2d , 

-=^«KK- 


( 11 ) 


This gives instead d = 12.55. 
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Adopting the range A - 16.63... 12.55 for predictions of the Witten-Sakai-Sugimoto model turns out to work 
remarkably well for a number of quantities of low-energy QCD at Nc - 3. While at infinite Nc, the anomaly of the 
Ua( 1) symmetry is suppressed and the rj' meson is a massless Goldstone boson, at finite Nc the model predicts a finite 
mass through a Witten-Veneziano formula (evaluated already in [34]) with the result 


niri' = 


1 

3 V37r 



( 12 ) 


With Mkk = 949 MeV and A ^ 16.63... 12.55 the numerical value for Nc — Nf — 3 turns out to be 967... 730 MeV, 
encompassing the experimental value 958 MeV (while the lower value also leaves room for a contribution due to finite 
quark masses that are not included in the chiral Witten-Sakai-Sugimoto model). 

Another quantity that is of particular interest with regard to glueball physics is the gluon condensate, which was 
calculated in Ref. [40] with the result 


c4 = (= 

7T ^ 


ANc 7 A 


(13) 


(see Ref. [25] for the precise translation to the conventions of ordinary QCD). For A = 16.63... 12.55 this yields 
= 0.0126.. .0.0072GeV^, which includes the standard SVZ sum rule value [41] at its upper end. 

In view of the aim of using the Witten-Sakai-Sugimoto model for a calculation of glueball decay rates, it is 
particularly interesting how its predictions for the decay rates of the p and a> meson compare with experiment. 

For the effective interaction Lagrangian of p and tt mesons one finds 


Zpnn = -gpnneabc(dpn‘‘)p^^‘n^, gpnn = 33.984 


(14) 


which yields 

Fp^ 2 ;rMp = 0.1535...0.2034 for 4= 16.63... 12.55, (15) 

including the experimental value [42] of 0.191(1). 

The decay co —> 3;7r, which requires the Chern-Simons part of the action (7), was calculated in Ref. [35]. The result 
is proportional to A~^ and therefore gives a larger span of predictions. 


F„^ 3 ;,/m„ = 0.0033...0.0102 for 4 = 16.63... 12.55. 


(16) 


The experimental value is [42] 0.0097(1). Encouragingly, this again fits the holographic result. 


GLUEBALL MASSES AND DECAY RATES IN THE 
WITTEN-SAKAI-SUGIMOTO MODEL 

In Ref. [36], the Witten-Sakai-Sugimoto models has for the first time been used for a quantitative evaluation of scalar 
glueball decay rates by including the glueball modes in the D8 brane action. In Ref. [38], we have repeated, partly 
corrected, and extended this work. 

With Mkk fixed to 949 MeV through the mass of the p meson, the predictions of the Witten model for the glueball 
spectrum are pinned down. The result is that the lowest scalar mode (5) has the mass Me — 855 MeV, which is only 
10% higher than nip and thus much lower than the predictions of quenched lattice simulations. However the lowest 
scalar mode in the tensor multiplet, which is a predominantly dilatonic mode, has the mass Mo - Mj - 1487 MeV, 
not far from lattice results as far as the O^'*' glueball is concerned, while the lattice predictions for the lowest tensor 
glueball is about 50% higher. 

The same picture arises if one considers the ratio of glueball mass over the square root of the string tension using 
Eq. (11) and compares with the lattice results of Ref. [24] obtained at large Nc- In the right panel of Eig. I, the 
holographic results are given by dotted green lines for the lowest scalar (exotic) mode, which is only half of the lattice 
result for the lowest O^'*' state, whereas the lowest dilatonic mode (dotted black) is slightly below the latter. (Note that 
“exotic” refers only to the polarization of the underlying graviton mode and not to an exotic assignment.) 

Although these results do not depend on the value of 4 within the Witten-Sakai-Sugimoto model, one can expect 
important corrections involving inverse powers of 4 and Nc in the form of higher-derivative terms from string- 
theoretic effects beyond the supergravity approximation. However, the discrepancy between Me and the lattice result 
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TABLE 1. Masses M and decay rates TQ^2n 
of lightest scalar and tensor glueballs 


Mode 

Mass 
M [MeV] 

TIM 

(4= 16.63... 12.55) 

0++ (E)* 

855 

0.092... 0.122 

0++ 

1487 

0.009... 0.012 

2++ 

1487 

0.015...0.019 


* “exotic” polarization 


is unusually large considering the other predictions of the model. As already remarked, the lowest scalar mode 
corresponds to an “exotic polarization” of the six-dimensional graviton in that it involves metric components SGaa 
referring to the extra spatial dimension whose only purpose in the Witten model is to implement the breaking of 
supersymmetry. Perhaps the exotic mode should be discarded from predictions of the Witten-Sakai-Sugimoto model, 
assuming that it will somehow disappear when the supergravity approximation is left in the limit Mkk —> oo, d —> 0 
that brings the model towards a string-theoretic dual of large-Ac QCD. 

In Ref. [36], decay rates have been studied only for the lowest (“exotic”) scalar mode. In Ref. [38], we have 
also worked out the interactions of the predominantly dilatonic mode, the tensor mode, and their excitations with 
pseudoscalar and vector mesons. Writing out here only the interactions with pions, the exotic scalar mode has an 
interaction Lagrangian of the form 


£G£- 


cidnTfdyTf —G 


Cl =62.655d“2 A:'M, 


-1 A/f-l 
KK’ 


Cl = 16.390 d“2A:‘M, 


-1 A/f-l 
KK’ 


(17) 


where we differ from Ref. [36] by a factor of V2 in the coefficient ci and by the term with the coefficient ci, which 
has been dropped in Ref. [36]. 

The predominantly dilatonic mode has an interaction with pions given by 




a _a.q _ci 
- Y dvTl 


rf^- 


d>^d' 


D, c/i = 17.226d“2 


-1 A/f-1 

KK’ 


( 18 ) 


and for the tensor mode we find 

^ ^1 ^ (i = Vhc/l, (19) 

where is normalized such that = \T^y{u-M^yT^" + + ..., with B,Bi^ being Lagrange 

multiplier fields. 

In Table 1 the resulting decay rates Tare given, which show that the relative width T/M is small, around 1 % and 
2%, for the dilatonic and the tensor mode, respectively. However, the result for the lowest (exotic) mode with mass 855 
MeV is an order of magnitude larger. (Parametrically, all three results are of the order AT^Ny^.) This peculiar result is 
in fact reminiscent of a scenario discussed by Narison and others in Ref. [46, 47, 48, 49], where both a broad glueball 
around 1 GeV and a heavier narrow glueball around 1.5 GeV were found to be required by QCD sum rules. There 
the lighter glueball was interpreted as a pure-glue component of the cr-meson /o(500). Recent lattice simulations of 
unquenched QCD [7] did not find changes of quenched result for the mass of the lowest scalar glueball around 1.5-1.8 
GeV, which we take as motivation to concentrate on the predominantly dilatonic mode and to investigate how its decay 
rates compare with experimental glueball candidates. 

In the range of 1.5-1.8 GeV, the Particle Data Group [42] lists two isoscalar O'*'^ mesons, /o(1500) and /o(1710). 
These are frequently and alternatingly discussed as possible manifestations of a large glueball component. In Ref. [38] 
we have evaluated the decay rate into four pions which are produced by glueball couplings with Tr{p{n,n'\) and Tr(pp). 
The dilatonic mode is also heavy enough to have decays into two kaons and two rj mesons. In the chiral Sakai-Sugimoto 
model, the latter are given simply by the flavor-symmetric factors 4/3 and 1 /3 compared to decay into two pions. 

In the upper half of Table 2, the resulting decay rates are compared with the experimental data for /o(1500) from 
Ref. [42]. We find that the decay into two pions comes out about 50% too low, while the decay into kaons is too large 
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TABLE 2. Experimental data for the decay rates of the 
isoscalar mesons /o(1500) and /o(1710) confronted with 
the holographic results for various decay channels of the 
predominantly dilatonic glueball (Go) in the chiral limit 
with glueball mass adjusted to the respective experimental 
values = 1505 and 1722 MeV, for ’t Hooft coupling 
varied from 16.63 to 12.55. 


decay 

r/M (exp.)* 

r/M[Gz)(M“P)] 

/o(1500) (total) 

0.072(5) 

0.027... 0.037 

/o(1500)^4;r 

0.036(3) 

0.003... 0.005 

/o(1500)^2;r 

0.025(2) 

0.009... 0.012 

/o(1500)^2A 

0.006(1) 

0.012...0.016 

/o(1500)^2;7 

0.004(1) 

0.003... 0.004 

/o(1710) (total) 

0.078(4) 

0.059... 0.076 

/o(1710)^2A 

I 0.041(2) 
^(0.047(17) 

0.012... 0.016 

/o(1710)^277 

10.020(10) 

*10.022(11) 

0.003... 0.004 

/o(1710)^2;r 

10.017(4) 

* 1 0.009(2) 

0.009... 0.012 

/o(1710)^4;r 

? 

0.024... 0.030 

/o(1710)^2^u^6;r 

seen 

0.011...0.014 


• Experimental data are from the Particle Data Group [42] except 
for those marked by a star (★), which are from Ref. [43] where the 
total width of /o(1710) was divided up under the assumption of a 
negligible branching ratio to four pions and two w mesons, using 
data from BES [44] (upper entry) and WA102 [45] (lower entry), 
respectively. 


by a factor of 2. The decay into four pions, which is the dominant decay mode of /o(1500), is missed by an order of 
magnitude. 

Turning to /o(1710), we need to manually adjust the mass of our dilatonic mode in order to take into account that 
this meson is above the threshold of 2p and 2a>. We do so by simultaneously rescaling the mass scale in the dilaton 
couplings, which leaves the dimensionless ratio T/M unchanged for the decays into two Goldstone bosons. However 
the decays into four pions and also two at mesons (which because of their narrow width are treated as nearly stable) 
now become significantly increased. Decay of /o(1710) into two co mesons has been seen according to the Particle 
Data Group [42], but not into four pions. The decay width into two pions [which is independent of the extrapolation to 
the higher mass of the /o(1710)] turns out to have the right magnitude, but the decay into two kaons is underestimated 
if only a flavor-symmetric enhancement of 4/3 is included. 

It is to be expected that a deformation of the Witten-Sakai-Sugimoto model towards finite quark masses and thus 
nonzero masses for the pseudo-Goldstone bosons such as discussed in Ref. [50] gives additional contributions to the 
vertices of a glueball with two pseudo-Goldstone bosons. This can significantly enhance the decays into two kaons and 
two T] mesons, as suggested by the lattice results of Ref. [51 ] and the analysis of Ref. [52], where this phenomenon was 
termed chiral suppression. The only feature of our results that would then deviate from the decay pattern of /o(1710) 
(as presently known) is the substantial decay rate of the holographic glueball Go into four pions. The experimental 
data for partial decay widths of /o(1710) given in Table 2 are from Ref. [43], where it was simply assumed that the 
partial widths for pairs of pseudoscalar mesons add up to the total width; in the analysis of Ref. [19], e.g., they add 
up to only 70%. Thus the experimental situation certainly needs further clarification, before one can judge if and by 
how much the holographic result overestimates the decay into four pions. If /o(1710) is a nearly unmixed glueball as 
suggested by Refs. [16, 18, 53, 54, 21], then the (extrapolated) prediction of the Witten-Sakai-Sugimoto model would 
be that there should a nonnegligible decay width into four pions and about half as much for the one into two w. (In 
Ref. [21], which employed an extended linear sigma model with a dilaton, the decay into four pions turned out to be 
strongly suppressed due to their very large value of the gluon condensate.) 

In Ref. [38] we have also evaluated the parametrically strongly suppressed decay of scalar glueballs into four tiP, 
which turns out to be significantly smaller than the (small) rate observed for/o(1500). This also supports the conclusion 
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that the result for the Witten-Sakai-Sugimoto model is not in agreement with a pure-glueball interpretation of /o(1500), 
but seems to prefer such an interpretation for /o(1710). 

For the tensor glueball we have obtained a very small decay width into two pions. But since the lattice points to 
a tensor glueball far above the threshold of two p mesons, decays into two vector mesons (possibly also including 
the (p) have to be considered as well. In Ref. [38] we have found that these decay channels dominate and lead to a 
comparatively large total width, larger than the /2 mesons listed by the Particle Data Group around and above 2 GeV, 
with the exception of /2(1950), which is indeed occasionally discussed as a glueball candidate. 


CONCLUSIONS AND OUTLOOK 

We have revisited the spectrum and the decay rates of glueballs in the Witten-Sakai-Sugimoto model which were 
first evaluated quantitatively in Ref. [36]. At variance with the latter we have concluded that only the predominantly 
dilatonic scalar glueball mode of the Witten model has a mass and decay rate comparable with glueball candidates in 
the range indicated by lattice gauge theory, while the lowest mode either has to be discarded or, more speculatively, 
could perhaps be identified with a broader glueball component of the cr-meson (cp. Ref. [46]) that only drops out in 
the theory without quarks. 

The decay pattern of the predominantly dilatonic scalar glueball is found to deviate rather markedly from the one 
observed for the glueball candidate /o(1500)—its dominant decay mode into four pions is not reproduced by the 
holographic results, while decay into two pions is underestimated by a factor 2. A somewhat better match is obtained 
for /o(1710), if this is assumed to be a nearly unmixed glueball [16, 18, 53, 54, 21], as the decay rate into two pions 
agrees rather closely with experimental data. The dominant decay into two kaons observed in experiments could be 
due to an enhancement of glueball coupling through mass terms of pseudo-Goldstone bosons, which will necessarily 
contribute to the coupling of glueballs to pseudoscalar mesons; we intend to study the effects of deformations of 
the Witten-Sakai-Sugimoto model away from the chiral limit in future work. Decay into four pions, however, are 
uncomfortably large when the mass of the holographic glueball is matched to that of /o(1710), given that this decay 
mode has not been observed in experiment. Decay into two a>, which mostly decay into three pions, has been observed. 
Our holographic result is that it should be at the level of half of the rate into four pions, and both not strongly 
suppressed. 

At the level our analysis, no mixing of glueballs with qq states is taken into account so that the results at best 
pertain to approximately unmixed glueballs. Mixing, which can strongly obscure the decay pattern of a pure glueball, 
is suppressed in the Witten-Sakai-Sugimoto model by at least a factor N~^ and would require further string-theoretic 
input to be studied semi-quantitatively. Absent that, it might be interesting to incorporate our holographic results in 
phenomenological models of mixing between glueballs and qq states. 
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